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ABSTRACT 

Context. Recent studies carried out with SOHO and Hinode high-resolution spectrometers have shown that the plasma in the off-disk 
solar corona is close to isothermal. If confirmed, these findings may have significant consequences for theoretical models of coronal 
heating. However, these studies have been carried out with diagnostic techniques whose ability to reconstruct the plasma distribution 
with temperature has not been thoroughly tested. 

Aims. In this paper, we carry out tests on the Monte Carlo Markov Chain (MCMC) technique with the aim of determining: 1) its 
ability to retrieve isothermal plasmas from a set of spectral line intensities, with and without random noise; 2) to what extent can it 
discriminate between an isothermal solution and a narrow multithermal distribution; and 3) how well it can detect multiple isothermal 
components along the line of sight. We also test the effects of 4) atomic data uncertainties on the results, and 5) the number of ions 
whose lines are available for the DEM reconstruction. 

Methods. We first use the CHIANTl database to calculate synthetic spectra from different thermal distributions: single isothermal 
plasmas, multithermal plasmas made of multiple isothermal components, and multithermal plasmas with a Gaussian DEM distribution 
with variable width. We then apply the MCMC technique on each of these synthetic spectra, so that the ability of the MCMC technique 
at reconstructing the original thermal distribution can be evaluated. Next, we add a random noise to the synthetic spectra, and repeat 
the exercise, in order to determine the effects of random errors on the results. We also we repeat the exercise using a different set of 
atomic data from those used to calculate synthetic line intensities, to understand the robustness of the results against atomic physics 
uncertainties. The size of the temperature bin of the MCMC reconstruction is varied in all cases, in order to determine the optimal 
width. 

Results. We find that the MCMC technique is unable to retrieve isothermal plasmas to better than AlogT ^ 0.05. Also, the DEM 
curves obtained using lines calculated with an isothermal plasma and with a Gaussian distribution with FWHM of log T ^ 0.05 are 
very similar Two near-isothermal components can be resolved if their temperature separation is A log T = 0.2 or larger Thus, DEM 
diagnostics has an intrinsic resolving power of log T = 0.05. Atomic data uncertainties may significantly affect both temperature and 
peak DEM values, but do not alter our conclusions. The availability of small sets of lines also does not worsen the performance of the 
MCMC technique, provided these lines are formed in a wide temperature range. 

Conclusions. Our analysis shows the present limitations in our ability to identify the presence of strictly isothermal plasmas in stellar 
and solar coronal spectra. 

Key words. Methods: data analysis — Techniques: spectroscopic — Sun: corona — Sun: UV radiation 



1. Introduction 

Detailed measurements of the temperature distribution of the 
plasma of solar and stellar coronae are of critical importance 
to understand how they are heated to multimillion degree tem- 
peratures, so it is no surprise that the thermal structure of the 
solar corona has been studied ever since the publication of the 
seminal paper of Edlen (1942), who discovered its high temper- 
ature for the first time. These studies have known a revival in 
recent years, thanks to the availability of a large body of obser- 
vations obtained with the new instrumentation carried on board 
the SOHO, TRACE, STEREO, Hinode and SDO satelUtes. 

Plasma loops are a fundamental component of the solar 
corona, since they are ubiquitous both in quiet and active so- 



lar plasmas. Steady-state models of loops predict that the tem- 
perature of these structures varies along the loop (Rosner et 
al. 1978), and that loops with different length can have very 
different temperatures; nanoflare-based models predict finely- 
structured loops filled with multithermal plasma all along their 
length (Patsourakos & Klimchuk 2005). A consequence of these 
models is that quiet coronal plasmas observed outside the so- 
lar limb should be multithermal because they are the sum of a 
multitude of loops with different lengths and temperatures along 
the line of sight. The question of clearly distinguishing between 
multithermal and isothermal plasma along the line of sight is 
thus very important in the framework of coronal heating models 
(e.g. Klimchuk 2006, Reale 2010). However, recent measure- 
ments have provided increasing evidence of a nearly isother- 
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mal nature of both coronal hole and quiet Sun plasmas when 
they are observed above the solar limb, contrary to model pre- 
dictions. Several authors (e.g. Feldman et al. 1998 and 1999, 
Warren 1999, Doschek et al. 2001, Landi 2008) reported single- 
temperature determination which are remarkably similar even 
if taken with different instruments and at different times dur- 
ing the solar cycle. A constant temperature was also found to 
characterize an entire streamer by Landi et al. (2006). Even a 
complex active region, when observed at the solar west limb, 
seemed to be composed by three near isothermal plasmas (Landi 
& Feldman 2008). Further studies carried out with Hinode/EIS 
observations, on the contrary, provided evidence that the quies- 
cent solar corona is not strictly isothermal, but it is characterized 
by a rather narrow temperature structure with a tail that extends 
to higher temperatures (Warren & Brooks 2009, Brooks et al. 
2009, Hahn et al. 201 1). The aim of the present paper is to as- 
sess how reliable such isothermal plasma claims are. 

These measurements have been caiTied out either with 
the Emission Measure (EM) loci technique or with tradi- 
tional Differential Emission Measured (DEM) techniques (see 
Feldman & Landi 2008 and references therein). The former al- 
lows to determine the temperature and the total emission mea- 
sure (EM) of an isothermal plasma, and it has been the tech- 
nique used in most of the studies that resulted in the isother- 
mal corona results. This technique, while useful in the case of 
strictly isothermal plasmas, fails at providing an assessment of 
the uncertainty of its results. Traditional DEM techniques rely 
on several different methods to reconstruct the thermal structure 
of multithermal plasmas determining the DEM curve. However, 
most of them rely on assumptions on the smoothness of the dis- 
tribution, that make it very difficult to study isothermal or near- 
isothermal plasmas. 

Landi & Klimchuk (2010) carried out tests aimed at evaluat- 
ing the ability of the EM loci technique to discriminate between 
isothermal and multithermal plasmas, when the spectral line in- 
tensities used for the reconstruction included uncertainties. They 
also developed a quantitative method to determine how multi- 
thermal could a plasma distribution be in order to be compatible 
with a set of observed line intensities that was also compatible 
with an isothermal condition. 

In the present work we extend the Landi & Klimchuk (2010) 
study to the application of DEM diagnostic techniques to spec- 
tral line intensities. In particular, we aim at: 1) determining the 
ability of such techniques to retrieve isothermal plasma from a 
set of spectral line intensities, with and without random noise; 
2) understanding to what extent are DEM diagnostic techniques 
able to discriminate between a truly isothermal plasma and a 
multithermal plasma with a naiTow distribution; 3) quantifying 
how well can such technique detect and resolve multiple isother- 
mal components along the lines of sight; 4) investigate the effects 
of atomic data uncertainties on the results; and 5) study how 
DEM reconstructions are affected by the number of ions whose 
lines are available. 

Similar studies with different goals were recently carried out 
in preparation for the SDO and Hinode missions by Golub et 
al. (2004) and Weber et al. (2004). Both studies were focused 
on determining how well, and with what uncertainties, could a 
DEM curve known a-priori be reconstructed using data from 
the limited number of broad-band filter intensities provided by 
Hinode/XRT and narrow-band filters in SDO/AIA. They also 
discussed how efficiently could DEM curves from the large num- 
ber of pixels expected from such images be determined and vi- 
sualized. Unlike the present work, such studies were aimed at 
imaging instruments and did not use spectral lines intensities; 



also, no guidelines were provided with regard to the interpreta- 
tion of almost isothermal results. Further, the results of such tests 
strongly depend on the specific temperature response function of 
each instrument, so that they can not be easily extended to other 
instruments. 

Line intensities can potentially provide better constraints on 
DEM curves due 1 ) to the better temperature resolution provided 
by their emissivities; and 2) to the finer sampling of the tem- 
perature range of the solar corona provided by the large num- 
ber of ions included in the spectral range of the available spec- 
trometers (CDS and SUMER on SOHO, EIS on Hinode). In 
the present work, we focus on spectral lines; given the larger 
amount of available constraints, we can hopefully obtain further 
insight on the intrinsic abilities of DEM diagnostic techniques at 
studying isothermal or near-isothermal plasmas. Our results can 
be applied to analyses of spectra observed by any of the avail- 
able spectrometers observing the solar upper atmosphere, such 
as Hinode/EIS, SOHO/CDS, SOHO/SUMER, SOHO/UVCS, 
FUNIS, SERTS, HRTS and, to a certain extent, even the lower- 
resolution SDO/EVE. In fact, different spectrometers working 
at different wavelength ranges may still observe lines from the 
same element, so that they will be sampling the same tempera- 
ture interval regardless of their passband. Thus, our results will 
be of more general use than Golub et al. (2004) and Weber et al. 
(2004). 

Many different DEM diagnostic techniques have been de- 
veloped in the past (see the reviews of Phillips et al. 2008 and 
Harrison & Thompson 1992). Harrison & Thompson (1992) car- 
ried out a comparative analysis aimed at determining which, 
among six different DEM diagnostic techniques, was most 
successful at reproducing pre-defined DEM curves using a 
set of emission lines. More recently, Mark Weber and Paul 
Boerner led a similar study in preparation to the SDO mission 
{http://www.lmsal.com/~boemer/demtest/), where, however, the 
results are not clearly discussed and no clear conclusions are 
drawn; also, such a study does not appear to have been pub- 
lished in the literature. Most importantly, Boerner & Weber did 
not address the ability at discriminating between isothermal and 
non-isothermal plasmas, which is the focus of the present work. 

After inspecting the material in the Boerner & Weber web- 
site, we found that Monte Carlo Mar-kov Chain (MCMC) method 
developed by Kashyap & Drake (1998) provided the most robust 
and accurate results. Also, this method dispenses from many 
of the assumptions common to other DEM diagnostic tech- 
niques, and provides an assessment of the uncertainties of the 
final DEM curve. Thus, we chose to focus on this method in 
the present study. Kashyap & Drake (1998) and Kashyap et al. 
(2004) briefly described tests made to ensure and assess the ro- 
bustness and reliability of the MCMC technique, but they did 
not focus on the ability of the MCMC technique to discriminate 
between an isothermal and multithermal plasma. Here we also 
provide a set of empirical guidelines for the choice of the input 
parameters and for the interpretation of MCMC reconstructions. 

Section |2] introduces the main plasma structure diagnostic 
techniques, including MCMC, and discusses their assumptions; 
it also explains the methodology we follow in our tests. Section[3] 
describes the results of our tests for several different ad-hoc dis- 
tributions, and Section|4]reviews the results of this study. 
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2. Analysis method 

2.1. Emission line intensities 

Solar coronal plasmas are usually tenuous enough to be optically 
thin below 2000 A. The intensity of an optically thin emission 
line can be written as 



X 



G(T,Ne)N,dV 



(1) 



where A^^ is the electron density, V is the emitting volume along 
the line of sight, d is the distance between the emitting source 
and the observer, and G{T,Ne) is the Contribution Function of 
the emitting line defined as 



G(T,Ne) 



where 



Nj(X^"') NiX^"') NiX) N(H) A 
N{X+"') NiX) NiH) Ne We 



(2) 



- is the relative population of the upper level j and 
depends on electron temperature and density; 

- ^^'^^^ ^ is the relative abundance of the ion X^"' (ion fraction); 
under ionization equilibrium conditions, it depends on the 
electron temperature; 

- is the abundance of the element X relative to hydrogen; 

- is the hydrogen abundance relative to the electron den- 
sity (=s 0.83 for fully ionized plasmas); 

- A is the Einstein coefficient for spontaneous emission. 

When the plasma is multi thermal Equation[T]can be rewritten by 
defining the Differential Emission Measure (DEM), ip(T), as 



G(T,Ne)^(T)dT with ^(T)^Nj^ 

di 



(3) 



The DEM indicates the amount of material in the plasma as a 
function of temperature. When the plasma is isothermal at tem- 
perature Tc, we can define the Emission Measure (EM) of the 
plasma as 



1 



G{T„Ne)EM with EM ■ 



r N^dv 

Jv 



(4) 



The EM of the plasma is a measure of the total amount of plasma 
at the temperature Tc. 

2.2. Tliermal structure diagnostic teciiniques 

When the plasma is isothermal at the temperature the EM 
can be determined from line intensities as: 



EM = And"- 



I 



G{Tc,Ne) 



(5) 



The EM loci diagnostic technique allows to measure simultane- 
ously the plasma EM and Tc values. It consists of calculating for 
each line the function EM(T) defined as 



EM{T) = And' 



G{T,Nc) 



EM(Tc) = EM 



(6) 



as a function of electron temperature, using the observed line in- 
tensities /. Since the plasma EM is the same for all lines, all 
the EM(T) curves, when displayed as a function of T in the 
same plot, should cross the same point coordinates (Tc,EM). 
The presence of this crossing point also confirms that the plasma 
is isothermal. When the plasma is not isothermal, the EM{T) 
curves do not cross the same point and the EM loci technique 
can not be used. 

When the plasma is multithermal, the plasma DEM needs 
to be determined. There are many different types of DEM di- 
agnostic techniques, and they have been reviewed by Phillips et 
al. (2008). Some techniques rely on the inversion of Equation[3] 
others on an iterative procedure; Monte Carlo methods are also 
available. 

The inversion of Equation[3]is carried out by discretising the 
temperature interval in bins whose width is chosen by the user. In 
each bin the DEM is assumed to be constant so that the observed 
intensity /, of the spectral line s can be written as 



1 

And^ 



1=1 



GAT)dT 



(7) 



The values of (p, are determined minimizing via an iterative tech- 
nique the sum the quantity H = S + x^, where S is the en- 
tropy of the the quantity tpi and has the standard definition. 
In this procedure, each observed line is arbitrarily associated to 
a temperature value Tg/f (discussed below), and the results for 
all lines whose Tg/f falls in the same bin are combined into a 
single value. The DEM curve is then determined by fitting the 
ipi values of all temperature bins with a polynomial or a spline 
function. 

Iterative techniques start from an arbitrary initial DEM curve 
and use each spectral line to calculate a correction C^Teff) asso- 
ciated to a temperature Tc/f (discussed below). c{Teffj values 
whose Teff lies in the same temperature bin (whose arbitrary 
width is chosen by the user) are averaged together. The result- 
ing average corrections for all bins are then interpolated with a 
polynomial or a spline function and applied to the initial DEM 
curve to calculate the new curve to be used in the next iteration 
as initial DEM. 

The definition of Tg/f for each line is rather arbitrary. In 
many cases Tc// is assumed to be T„iax, the temperature of max- 
imum abundance of the ion emitting the line. In other cases, 
it is defined as some sort of DEM- and Contribution Function- 
weighted mean of the temperature, such as in Landi & Landini 
(1997). 

2.3. Pitfalls of the diagnostic techniques 

Both the inversion and iterative DEM diagnostic methods rely 
on three main assumptions: 1) the interpolation with a polyno- 
mial or a spline function assumes implicitly that the plasma is 
multithermal and may oversmooth the results; 2) the measured 
DEM depends on the widths of the temperature bins chosen to 
group the lines; and 3) with only few exceptions, Tcjj is arbi- 
trarily chosen as T^ax, regardless of the thermal structure of the 
plasma. 

The first problem makes the DEM analysis unsuitable for 
determining whether a plasma is isothermal or multithermal, be- 
cause the DEM analysis implicitly assumes the latter condition. 
Also, associating each spectral line to a predetermined temper- 
ature like Tmax is likely to bias the results against an isothermal 
solution: in fact, if the plasma has a very narrow temperature 
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distribution, the temperature of the plasma may be very different 
from the r,„,„ of most Hnes in the data set. Associating each Hne 
to Tmax results in associating a line to a temperature where there 
might not be any plasma at all and this leads to an artificially 
multi thermal DEM. 

If the width of the temperature bin is chosen too large, vari- 
ations of the DEM with temperature occurring in temperature 
intervals smaller than the bin width are averaged out so that the 
final DEM is oversmoothed. If on the contrary the bin width is 
too narrow, spurious effects due to problems in individual lines 
and photon noise can severely affect the solution. 

Moreover, the EM and DEM techniques require the use of 
lines emitted by a large number of different ions whose ion frac- 
tions are non-negligible in widely different temperature ranges. 
The use of only few lines greatly limits the diagnostic capabil- 
ities of these techniques, as demonstrated by Judge (2010) and 
Landi & Khmchuk (2010). 

2.4. The Markov Chain Monte Carlo DEM diagnostic 
technique 

The problems associated with DEM interpolation and with the 
definition of T^ff do not affect the Markov-Chain Monte Carlo 
(MCMC) technique developed by Kashyap & Drake (1998). 
This technique is based on a Bayesian statistical formalism that 
allows the determination of the most probable DEM curve that 
reproduces the observed line intensities. The heart of this tech- 
nique relies on the application of the Bayes theorem to line in- 
tensities, stating that the probability P(X, F) of obtaining a set of 
observed line intensities F = {F\,F2, F„) from a DEM char- 
acterized by a set of parameters X - {X\ , X2, X„) is given by 

P(X, F) = P(X) (8) 

where P{F) is a normalization factor, P(X) is an a priori proba- 
bility of the set of parameters X, and P(X, F/) is the probability 
of obtaining the observed intensity F, with the set of parameters 
X and has been defined by Kashyap & Drake (1998) as 

P(X,Fd^e-i(''-''"')'^''^^f (9) 

where Ff is the intensity of line / calculated with the DEM de- 
scribed by the set of parameters X, and cr, is the uncertainty of 
the observed intensity Fj. This method is implemented by choos- 
ing a grid of temperature bins and assuming that within each 
bin the plasma is isothermal and can be described by an emis- 
sion measure value EM/. Once the bin size is chosen, the set of 
parameters that describe the DEM is thus X = {EMi , EMm). 

In order to determine the set of parameters X that provide 
the maximum probability P(X,F), Kashyap & Drake (1998) 
adopted a Markov-Chain Monte Carlo approach. In this ap- 
proach, the set of parameters X that describes an initial trial 
DEM is varied step by step; in each step only one of the parame- 
ters in the set X (i.e. one EM, value only) is varied and the others 
are left unchanged; the change introduced to the varied parame- 
ter only depends on the set of parameters X of the previous step. 
The new set of parameters X' has a different probability P(X', F) 
from the previous one, and is accepted or rejected according to 
the Metropolis algorithm (Metropolis et al. 1953), based on the 
change of probability P. This algorithm consists of generating a 
random number u, such that < u < 1, and a function A(X,X') 
defined as 

[ P(X\F) 

A(X,X') = min 1, — — ^ (10) 
V ' ; P(X,F) 



The new set of parameters is accepted if m < A(X,X'), and 
rejected otherwise. In this way, not only is a new set of pa- 
rameters X' with greater probability than the previous one al- 
ways accepted (since A(X,X') - 1), but also a new set with a 
smaller probability has a finite chance of being accepted. This 
latter property helps finding the best distribution of parameters 
X by moving the solution out of local maxima. When the system 
has found the best solution, the distribution of the EM values in 
each bin found in all steps is used to determine the confidence 
interval at each temperature, thus providing an estimate of the 
uncertainty of the DEM at all temperatures. 

This method dispenses from two of the assumptions com- 
mon to the other two methods. First, it is not affected by the 
choice of T^jj, since each line contributes to determine the so- 
lution in the entire temperature range where its G(T) function is 
defined. Also, no interpolation of the solution is required and no 
assumption on the smoothness of the solution is made. However, 
the choice of the bin width is crucial for this technique as it is 
for the iterative and inversion ones, and it is necessary to assess 
its effect. 

2.5. Methodology of the tests 

In order to test the ability of the MCMC method to recon- 
struct the DEM, we apply this technique to line intensities calcu- 
lated using known thermal distributions and version 6.0. 1 of the 
CHIANTI database (Dere et al. 1997, 2009) using the Bryans 
et al. (2009) ion fractions; thus, uncertainties due to the atomic 
physics do not affect our tests. We assumed a density value of 
\ogNe - 9.0 (Ne in cm"^), typical of moderately active region 
plasmas. In order to understand how the MCMC diagnostic tech- 
nique is affected by errors in the data (both in the atomic data 
used for emissivity calculations as well as observational uncer- 
tainties), for each of the synthetic spectra we used we generated 
five additional datasets; in each additional dataset we have added 
random errors of up to 20%, i.e. we have randomized each line 
intensity within a 20% range around the nominal value. Thus, 
each thermal distribution we have considered will be recon- 
structed using first the original synthetic line intensities with no 
errors added, and then each of the five different sets of intensities 
with random errors. With five different "noisy" datasets, the sta- 
bility of the solution against errors in the data can be evaluated. 
As a further test on the effects due to atomic physics only, we re- 
peated the MCMC determinations using the synthetic intensities 
calculated with version 2 of the CHIANTI database (Landi et al. 
1999) and the ion fractions of Mazzotta et al. (1998). The atomic 
data and ion fractions in the two CHIANTI versions are very dif- 
ferent, so that the use of CHIANTI V.2 represents a rather bad 
case of uncertainties in atomic parameters. 

We considered the particularly favorable situation of many 
lines available for diagnostics: we used lines from 45 differ- 
ent ions, listed in Table [T] that allow to sample very finely the 
temperature interval 4.6 < logT < 8.0, as shown by Figure [T] 
However, above log T = 7.0 and below log T = 5.5 the number 
of available ions decreases, and thus the MCMC reconstruction 
is expected to be less accurate. 

The arbitrary width of the temperature bin has been varied, 
in order to determine the optimal value that allows us to best 
reconstruct the initial thermal distribution and its temperature 
dependence while limiting both the noise and variability due to 
a too small bin width and the oversmoothing due to a too large 
width. We run our tests using a grid of width W values: W = 
0.01,0.02,0.05,0.1,0.2. 

We considered several types of thermal structure: 
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Ion 


Wvl.(A) 


log r„„_v 






Ion 


Wvl.(A) 


l0§ Tmax 


G(r,„nv) 




On: 


702.900 


4.90 


1.13e-13 




Six 


871.726 


6.05 


1.97e-16 




Oiv 


787.710 


5.17 


l.OOe-12 




Sx 


264.231 


6.18 


7.31e-15 


★ 


Ov 


629.732 


5.37 


7.60e-12 




Sxi 


281.402 


6.28 


2.45e-15 


★ 


Ovi 


1031.914 


5.48 


2.91e-12 


★ 


Sxii 


288.421 


6.35 


4.76e-15 


★ 


Neiv 


543.887 


5.21 


6.53e-14 




S XIII 


256.685 


6.42 


1.75e-14 


★ 


Ne V 


572.336 


5.42 


l.OOe-13 




Fe vm 


185.213 


5.62 


9.55e-14 


★ 


Ne VI 


558.685 


5.59 


1.02e-13 




Feix 


171.073 


5.87 


5.24e-13 


* 


Ne VII 


465.220 


5.71 


5.47e-13 




Fex 


177.240 


6.04 


1.12e-13 


* 


Ne VIII 


770.410 


5.80 


2.54e-13 




Fexi 


182.169 


6.13 


2.87e-14 


★ 


Mg V 


276.579 


5.45 


2.51e-14 


★ 


Fexii 


193.509 


6.19 


9.84e-14 


★ 


]VIg VI 


Z lyt.Jyl 


J.OJ 


z. 1 le- iH- 


"At 


Fe XIII 






0. Joe- m- 


"At 


Mg vii 


276.154 


5.78 


9.33e-15 


★ 


Fexiv 


264.790 


6.29 


5.79e-14 


* 


Mg VIII 


782.364 


5.90 


6.18e-15 




Fexv 


284.163 


6.34 


3.18e-13 


* 


Mgix 


706.036 


5.99 


1.95e-14 




Fexvi 


262.976 


6.43 


l.lOe-14 


★ 


Mgx 


609.794 


6.07 


1.54e-13 




Fe XVII 


204.665 


6.61 


9.20e-16 


* 


Si VI 


246.003 


5.61 


2.91e-14 


★ 


Fe XVIII 


974.860 


6.86 


7.91e-15 




Si VII 


275.361 


5.78 


6.15e-14 


★ 


Fexix 


1118.057 


6.95 


5.14e-15 




Si VIII 


276.850 


5.93 


1.13e-14 




Fexx 


132.840 


7.01 


7.55e-15 




Si IX 


258.082 


6.05 


3.18e-15 


★ 


Fexxi 


1354.067 


7.06 


1.27e-14 




Six 


258.371 


6.15 


3.78e-14 




Fe XXII 


845.571 


7.11 


6.10e-15 




Si XI 


580.920 


6.22 


9.08e-15 




Fe XXIII 


263.766 


7.17 


1.73e-15 




Si XII 


520.666 


6.30 


4.48e-14 




Fe XXIV 


192.029 


7.26 


1.17e-14 


★ 


Sviii 


198.554 


5.91 


3.93e-15 















Notes. r„,„,.: temperature of maximum ion abundance. G(T„ax)'- hne emissivity at r„,„,.. Stars indicate the ions whose lines can be observed by the 
Hinode/EIS spectrometer (although EIS lines may be different from those listed here). 



3. A single Gaussian distribution of variable width; 

These thermal distributions are aimed at determining 1) to what 
extent and with what precision is the MCMC diagnostic tech- 
nique able to detect an isothermal plasma; 2) how capable is it 
at separating two isothermal components close in temperature; 
3) to what extent is it able to discriminate between an isother- 
mal component and a multithermal plasma characterized by a 
Gaussian distribution; and 4) how do random errors and atomic 
physics uncertainties affect all the above results. 



2.6. Uncertainties of fhe results 

One of the main advantages of the MCMC technique is the abil- 
ity to provide an estimate of the confidence level of the solu- 
tion. To determine it, we used the following procedure. Once the 
best solution is found, we have run additional 50,000 runs of the 
Markov Chain procedure allowing each bin to randomly change 
by factors up to 6 orders of magnitude. For each change, we 
calculated the of the solution and its probability, defined by 
Equation|9l We arbitrarily retained only the modified EM distri- 
butions with a probability of at least 0. 1 xPtesi, where Pbest is the 
probability of the best solution. Then, for each temperature bin 
we have calculated the average log EM value and its variance, 
and taken the value of the latter as an indication of the uncer- 
tainty on the EM value of that bin in the best solution. 

2.7. Evaluation of the results 

The resulting EM(r) curve is compared to the EM distribu- 
tion used to calculate the synthetic spectra. The parameters we 
used to assess the quality of a reconstruction are shown in 




Fig. 1. Temperature coverage of the ions listed in Table [T] For 
each ion only the temperature interval where the fractional abun- 
dance is larger than 0.01 has been considered. Ion abundances 
are from Bryans et al. (2009). Top: number of ions emitting in 
each temperature bin (bin width W = 0. 1). Bottom: Temperature 
range covered by the ions of each element. 



1. Single isothermal plasma; 

2. Two isothermal plasmas, with variable peak temperature sep- 
aration and variable relative peak amplitude; 
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Table |2] In this Table, we report the cases where the plasma 
is isothermal with log T - 6.0, or made as the sum of two 
isothermal components with temperatures logTi - 6.0 and 
log T2 - log T\ + A log T . The temperature distance of the two 
isothermal components in the multithermal case is chosen to be 
Alogr = 0.10,0.15,0.20. A successful fit to the original DEM 
curve is evaluated using the following criteria: 



1. Sms\\x^\ 

2. the DEM(T) curve should have a single peak (or two peaks 
for the multithermal case) at the temperatures log T\ and 
log 7^2 ; 

3. there should be no spurious components far from T\ and T2; 

4. the peak should be unique, i.e. there should not be spurious 
peaks anywhere near Ti and T2 (noise in the solution); 

5. the width of the peak (taken as the temperature range around 
Ti and T2 where the EM(T) value is 1/10 of the peak value 
or larger) should be as small as the temperature bin W. 



3. Results 

3.1. Isothermal plasmas 

The isothermal plasma was chosen to have log T = 6.0 (in K) 
and log EM = 43.0 (in cm"-'). This temperature was chosen 
as it allowed us to include lines routinely observed by avail- 
able specti-ometers such as SOHO/CDS, SOHO/SUMER, and 
Hinode/EIS. Figure|2]shows the EM loci technique applied to the 
calculated intensities, without random errors included: the origi- 
nal temperature and EM values are recovered without problems. 
Figure[3]shows the MCMC reconstruction obtained without ran- 
dom errors and varying the bin size W. Results are summarized 
in Table |2l In all cases the MCMC technique is able to provide 
an isolated peak several order of magnitude larger than the back- 
ground. In the case of bin width W - 0.01 the shape of the peak 
is irregular and one additional satellite peak is present at higher 
temperatures: this feature is an artifact and it disappears when 
the bin width increases. In all cases the peak maintains a non- 
negligible width A(log T) X! 0.03 - 0.05 that is approximately 
constant as the bin width W increases. We interpret this as the 
ultimate resolution capability of the MCMC technique to resolve 
a purely isothermal plasma: the result is indistinguishable from a 
multithermal plasma with a very narrow width. Even if the;^'^ of 
the solution with W=0. 1 is the lowest, the solution with W=0.05 
is to be prefeiTed because the bin width is closer to the appar- 
ently intrinsic resolving power of the MCMC technique. In the 
case with W = 0.05 the two bins at the side of the peak show 
EM values of » 1 % of the peak value, and might be taken as real. 
Therefore, care should be taken in interpreting such features. 

The solutions provided by the MCMC technique always 
have a background at very low EM values at all temperatures 
(log iiM ^ 33 - 37). This background represents the maxi- 
mum values of the EM at all temperatures outside the isother- 
mal peak that provide contributions to the smaller than the 
computer precision. However, this background has no physical 
significance. 

Figure |4] shows how the results change when the random er- 
rors are added to the calculated intensities. The plasma EM ob- 
tained with W = 0.05 is displayed as a full black line, and the 
EM curves obtained by the five different "noisy" datasets are 
shown in red. The main peak remains unaltered, and since the 
wings, always present, are much lower than the peak value the 
technique is still able to recover the single peak. However, the 
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Fig. 2. EM loci technique applied to the single isothermal 
plasma, without random errors. 
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Fig. 3. EM(T) reconstruction of a single isothermal plasma, 
without random errors. The width of the temperature bin changes 
from W = 0.01 to 0.10. 
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Fig. 4. EM(T) reconstruction of a single isothermal plasma, 
with random errors within 20% added to line intensities (see text 
for details). The width of the temperature bin is W = 0.05. All 
five "noisy" datasets are shown. 
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W=0.01 


W=0.02 


W=0.05 


W=0.1 




Isothermal plasma - 


no random error 




0.013 


0.037 


0.0018 


0.00063 


Spurious comp. 


Y 


N 


N 


N 


Single peak 


N 


N 


Y 


Y 


Noise 


Y 


N 


N 


N 


EM(T) width 


0.03 


0.04 


0.05 


0.10 


log T\ 


5.99-6.03 


5.98-6.02 


6.00 


6.00 




Multithermal plasma 


- no random error 






A log r = 


: 0.10 






0.0062 


0.0085 


0.0058 


0.0069 


Spurious comp. 


Y 


N 


N 


N 


Two peaks 


N 


N 


Y 


Y 


Noise 


Y 


Y 


N 


N 


EM(T) width 


0.03 


0.06 


0.05 


0.10 


log 7^1 


5.99-6.02 


5.98-6.04 


6.00 


6.00 


log 7^9 


6.10-6.11 


6.10 


6.10 


6.10 






AlogT = 


= 0.15 






0.0013 


0.0035 


0.00087 


0.23 


Spurious comp. 


Y 


Y 


N 


N 


Two peaks 


N 


N 


Y 


Y 


Noise 


Y 


Y 


N 


N 


EM(T) width 


0.03 


0.06 


0.05 


0.20 


log 7^1 


5.99-6.02 


5.98-6.04 


6.00 


6.00-6.20 


lo? 


6.15 


6.14-6.16 


6.15 


6.00-6.20 






Alogr = 


= 0.20 






0.0053 


0.0078 


0.0017 


0.00055 


Spurious comp. 


Y 


N 


N 


N 


Two peaks 


N 


N 


Y 


Y 


Noise 


Y 


Y 


N 


N 


EM(T) width 


0.04 


0.04 


0.05 


0.10 


logTi 


5.99-6.02 


5.98-6.02 


6.00 


6.00 


log 7^2 


6.18-6.22 


6.18-6.22 


6.20 


6.20 



Table 2. Performance of the MCMC technique when appUed to 
isothermal plasmas, and to multithermal plasmas composed of 
two different isothermal components. See Section lZTl for details. 



background noise now is more significant and secondary peaks, 
much lower than the one at log T - 6.0, are sometimes present 
such as those at log T = 6.30 and 6.45. 

3.2. Two-component plasmas 

Figure |5]displays the results of the reconstruction of the EM(T) 
curve in the case of two isothermal components with the same 
EM (log 43.0 each, in cm"^, shown as red diamonds) and with 
variable temperature separation A log T and no random errors. 
The bin width W in each reconstruction was kept smaller than 
the temperature separation, and when possible we experimented 
with increasing W to determine how it affected the ability of 
the technique to separate the two components. We also studied 
AlogT - 0.20 and 0.25, but the results are no different than 
those obtained with AlogT - 0.15. Table |2] reports the results 



DEM reconstruction of isothermal plasmas 
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Fig. 6. EM(T) reconstruction of a a plasma made of two isother- 
mal components whose peak temperatures are separated by 
D(log T) = 0.15 (top) and 0.20 (bottom). The bin width W is 
0.05. Line intensities have been added a random errors within 
20% (see text for details). The reconstruction has been made 
without (black full curve) and with (red dashed curves) a ran- 
dom errors within 20% to each line intensity. 



for the reconstructions obtained with AlogT = 0.10,0.15 and 
0.20. 

Figure |5] shows that even when random errors are absent the 
MCMC technique is unable to resolve the two components when 
their temperature separation is smaller than A log T - 0.10, and 
provides a single, very noisy peak with a larger width. Also, 
small values of the bin width W worsen the noise in the solu- 
tion, further preventing the separation of the two components. 
At A log T = 0. 10 a complete separation of the two peaks can be 
achieved only when W - 0.05, because noise confuses the result 
for smaller bin widths and the presence of high EM values in the 
only bin between the peaks prevents definitive conclusions about 
separate components. Only when Alog T = 0.15 the two com- 
ponents can be fully and convincingly resolved with W = 0.05. 
Noise in the definition of each peak is significant at any value of 
ALogT for W < 0.05; the bin W = 0.1 allows us to resolve the 
two peaks only when A log T = 0.20 and causes trouble if the 
separation is Alog T - 0.15, providing a large 

The presence of random eiTors further limits the ability of 
the MCMC technique to resolve the two components. Figure |6] 
shows the MCMC reconstruction for A log T = 0.15 and 0.20 in 
the presence of random errors, with W - 0.05. The two peaks 
are still visible but badly resolved when Alog T - 0.15, while 
they are well separated in the other case. Thus, we conclude that 
AlogT = 0.20 is a reliable estimate of the smallest temper- 
ature distance between two isothermal components which the 
MCMC technique can realistically resolve when the lines listed 
in Table [T] are available. If only a smaller number of lines can 
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Fig. 5. EM(T) reconstruction of a a plasma made of two isothermal components whose peak temperatures are separated by £>(log T), 
which is varied from 0.025 to 0. 15. The bin width W, reported in each panel, is chosen to be always smaller than the peak separation 
D(logr). 



be used, the minimum A log T value may be significantly larger. 
It is also important to note the spurious components that arise 
around the two isothermal components when the random errors 
are present (e.g. at log T - 5.7, 5.8, 6.35 etc). Their peak values 
are significantly lower than the true components, but neverthe- 
less they complicate the interpretation of the result. 

The MCMC technique maintains its ability to resolve dou- 
ble peaks even when their relative size is much different. Tests 
have shown that when the peak separation is A log T - 0.20 the 
two peaks can still be resolved in the presence of random er- 
rors even if one of them is a factor 10 or even 100 smaller than 
the other as shown in Figure [7] where MCMC reconstructions 
of a plasma distribution with two isothermal components sepa- 
rated by A log T = 0.20 with different heights are shown without 



(black full lines) and with (red dashed lines) the presence of ran- 
dom errors. The peaks can be recognized because the peak EM 
values are very close in all reconstructions, while in the other 
temperature bins the EM values are scattered over a broad EM 
range. 

3.3. Gaussian DEM distribution 

The reconstruction of the Gaussian EM(T) is shown in Figure[8] 
The Gaussian EM{T) curves were defined with a variable 
FWHM (in log r): 0.05, 0.10, 0.15. The EM{T) curves ai-e also 
displayed in Figure |8] starting from top (FWHM=0.05) to bot- 
tom (FWHM=0. 15). The reconstruction was performed with dif- 
ferent widths W = 0.01,0.02,0.05. When W < 0.05, the solu- 
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Fig. 7. EM{T) reconstruction of a plasma made of two isother- 
mal components whose peak temperatures are log T = 6.0 and 
6.2. The hotter peak EM is a factor 10 (top) and 100 (bottom) 
smaller than the colder one. The bin width W is 0.05. The recon- 
struction has been made without (black full curve) and with (red 
dashed curves) a random errors within 20% to each line inten- 
sity. 



tions are too noisy to be acceptable in all cases, and also provide 
a solution too narrow when FWHM=0.05. When W - 0.05, the 
reconstruction is smoother and follows more closely the original 
EM(T) curve. 

When FWHM=0.05, the solution provided by the MCMC 
technique is very similar to the one provided in the single 
isothermal plasma component. Figure |9] compares the isother- 
mal results for the isothermal plasma component (black) and the 
Gaussian EM(T) (red): the two reconstructed curves are almost 
identical for FWHM=0.05, while the Gaussian reconstruction is 
significantly larger than the isothermal one when FWHM=0. 10. 
This means that the MCMC curve has problems distinguishing 
an isothermal plasma from a Gaussian plasma when the width 
of the latter is very small. This result means that the MCMC 
diagnostic technique is unable to discriminate between isother- 
mal and narrow multithermal distributions. If we combine this 
limitation with the intrinsic width found in the single, isother- 
mal plasma results, we can state that the MCMC technique is 
unable to unambiguously detect an isothermal plasma, and that 
is at best able to determine that a plasma distribution can have a 
log T width of 0.05 or smaller. 

3.4. Effects of a smaller number of Ions 

The tests that we have carried out so far included lines from 
forty-five ions, much more than those usually available in 
real-case scenarios of high-resolution observations carried out 
with existing instruments like Hinode/EIS, SOHO/CDS, and 




Fig. 9. Comparison of MCMC reconstruction of a single isother- 
mal plasma (black) and a Gaussian EM(T) curve (red). The 
FWHM of the Gaussian is 0.05 (left column) and 0.10 (right 
column). The temperature bin W is 0.01 (top row), 0.02 (middle 
row) and 0.05 (bottom row). 



SOHO/SUMER. These tests show the intrinsic strengths and 
weaknesses of the MCMC technique. However, the number of 
spectral lines and the amplitude of the temperature range they 
are formed in are also limiting factors in the overall accuracy 
of DEM diagnostic techniques, as shown by Judge (2010) and 
Landi & Klimchuk (2010). The question is, how do the results 
change when a smaller number of ions is available. 

To investigate this, we considered the most recent high- 
resolution spectrometer available for studies of the solar atmo- 
sphere: Hinode/EIS (CuUiane et al. 2007). We carried out the 
same tests described in the previous sections using only lines 
from the ions that can be observed by EIS. These are marked by 
asterisks in Table [1] and they are a total of 24 ions; we have omit- 
ted Si vm and S viii since their EIS lines are too heavily blended 
with other species to be used for DEM analysis. On the over- 
all, the results that we have reached in the previous sections are 
confirmed. Often, even smaller groups of ions than those avail- 
able in the EIS wavelength range are used. This is due to the fact 
that when EIS observations are designed and carried out, instru- 
mental or observational constraints usually make it necessary to 
transmit to the ground only portions of the EIS spectral range, 
and to sacrifice several ions. A typical dataset is the one used by 
Brooks et al. (2009), where the ions used were Si vii. Six, S x, 
and Fe viii-xv, for a total of 1 1 ions. Similar numbers are typ- 
ical of most SOHO/CDS and SOHO/SUMER observations as 
well. Results obtained from such a restricted dataset do not have 
significant differences when dealing with plasmas made of one 
or more isothermal components, and the MCMC technique is 
able to provide the same results as with larger datasets, provided 
that the plasma temperature lies inside the range of formation of 
these ions. 



3.5. Effects of Inaccuracies In atomic data 

We test how systematic errors due to uncertainties in atomic 
physics affect the DEM diagnostics by creating an inconsistency 
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Fig. 8. MCMC reconstructions of Gaussian EM(T) distributions. Original curves are also shown. The temperature bin width W is 
0.01, 0.02 and 0.05 (left, middle and center column, respectively). The Gaussian FWHM is 0.05 (top row), 0.10 (second row), and 
0.15 (bottom row). 



between the generation and the reconstruction of the test DEM. 
In particular, we use an older version (V.2) of the CHIANTI 
atomic code to calculate the contribution functions to be used 
for the DEM reconstruction. 

Figures[TO]and[TT]display (in blue) the results obtained when 
the emissivities of CHIANTI version 2 are used to determine 
the EM{T) curves using lines calculated with CHIANTI version 
6.0.1. The red lines indicate the results obtained in the previ- 
ous sections, for comparison purposes. Figure [TOl illustrates the 
isothermal and single-Gaussian (with W - 0.05) cases: a single, 
almost identical peak is retrieved in both cases, confirming that 
the DEM method is still unable to distinguish between a true 
isothermal plasma and one with a narrow Gaussian distribution. 
However, one difference to be noted is that the peak of the single 
component is slightly shifted towards lower temperatures. This 
is mostly due to the Mazzotta et al. (1998) Fe ion abundances 
used with CHIANTI V.2 being shifted towards lower tempera- 
tures than the Bryans et al. (2009) ones used with CHIANTI 
V.6. The level of noise in the V.2 solution is comparable to that 
of the datasets with random noise added. 

Figure [TT] shows the case of two isothermal components, 
with separation AlogT - 0.15 and 0.20. In the latter case, the 
two peaks are fully retrieved and separated, and the level of noise 
of the solution is comparable to the noise in the solutions ob- 
tained when random errors are included in the simulated spectra. 
When AlogT - 0.15, the two peaks are less clearly separated. 
The temperatures of the two peaks are slightly lower than in the 
original model, but their separation is coiTect. The only marked 
difference between the solutions obtained with V.6 and V.2 emis- 
sivities lies in the value of the peak EM in the hot isothermal 



component, since V.2 emissivities underestimate this by a factor 
5 and 2 in the A log T - 0.15 and 0.20 cases, respectively. 

4. Conclusions 

In this work we have tested the MCMC diagnostic technique by 
Kashyap & Drake (1998) in order to: determine its capability 
to reconstruct a given plasma thermal distribution (isothermal 
or multicomponent) from a set of individual line intensities; to 
study how the presence of random errors affects the MCMC re- 
constructions; to investigate the optimal size of the temperature 
bin width W as a compromise between retaining as much tem- 
perature structure as possible and minimizing the effect of noise 
in the solution; and to test the effect of different atomic data sets 
on the final results. We did this by applying the MCMC tech- 
nique to sets of lines whose intensities were calculated using a 
known thermal distribution first, and then randomized the inten- 
sities within 20% to simulate experimental uncertainties. 

Our study shows that the optimal bin width W is W = 0.05, 
as smaller values cause the MCMC technique to provide spu- 
rious results and a high level of noise in the final reconstruc- 
tion, while larger values oversmooth the results. Even when the 
emitting plasma is strictly isothermal, the MCMC technique is 
unable to distinguish between a truly isothermal solution and 
a Gaussian DEM with FWHM=0.05; also, the MCMC tech- 
nique is able to separate multiple near-isothermal EM compo- 
nents only if their separation in temperature is A log T - 0.20. 
Atomic data uncertainties can affect the results by providing less 
accurate peak EM values and shifting the EM peak temperature, 
but the A log T - 0.20 resolving power of the MCMC technique 
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Fig. 10. Comparison of EM(T) reconstructions of a plasma 
made of a single isothermal component (top) and a single 
Gaussian component (bottom); each reconstructed using two dif- 
ferent versions of CHIANTI: version 6.0.1 (red curves) and 2 
(blue curves). Thin lines in the top panel indicate the case with 
random noise added, from Section lTTI 



Fig. 11. Comparison of EM(T) reconstractions of a plasma 
made of two isothermal components separated by AlogT = 0.15 
(top) and AlogT - 0.2 (bottom); each reconstructed using two 
different versions of CHIANTI: version 6.0.1 (red curves) and 2 
(blue curves). Thin lines in the top panel indicate the case with 
random noise added, from Section ITTI 



is unaffected. The number of available ions does not affect the 
quality of the reconstruction, provided these ions are formed 
over a temperature range larger than the range where the plasma 
EM is significant. A smaller ion formation temperature range de- 
creases the temperature resolution achieved by the MCMC tech- 
nique. 
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